clc,clear;
x=1900:0.01:2300;
y=10:0.0001:14;
for i=1:40001
Z{1}(i)=Theorem2(1871.2,0.53512,0.30758,1800,x(i));
Z{2}(i)=Theorem2(9.64,0.01643,0.163104,9,y(i));
end
Z{1}=Z{1}';
Z{2}=Z{2}';
% plot(x,Z{1},'r.')
% legend('\mu=0.53512,\sigma=0.30758')
% title('the density function of sup_{0<t<La}U(t) for Shanghai Composite Index')
% xlabel('sup_{[0,La]}U(t)'),ylabel('Probability Density');
plot(y,Z{2},'r.')
legend('\mu=0.01643,\sigma=0.163104')
title('the density function of sup_{0<t<La}U(t) for Greatoo Intelligent Equipment Inc.')
xlabel('sup_{[0,La]}U(t)'),ylabel('Probability Density');
